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1. Introduction 

One of the current challenges in theoretical particle physics is to compute properties 
of the strongly coupled sQGP discovered at RHIC. AdS/CFT tools have given us 
some insight into the strongly coupled thermodynamics of gauge theories [1, 2, 3, 4, 
5, 6, 7]. However, it remains a mystery why these, mostly M = 4 supersymmetric, 
YM calculations work well for QCD. Part of the challenge is to either understand 
why this is so, or to find AdS duals of theories resembling QCD closer than Af = 4 
SYM. In this latter context a phenomenological AdS dual to Chiral perturbation 
theory constructed by Erlich et.al. is perhaps a good candidate [8]. It predicts the 
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meson spectrum within 10%-15% accuracy; and with a smooth instead of a hard 
IR-cutoff predicts the critical T c surprisingly well [9, 10]. 1 

This IR-cutoff is the essential new ingredient in AdS / QCD compared to AdS / CFT. 
Here we shall investigate the effects of this cutoff on photon and dilepton produc- 
tion rates at strongly coupling. Remarkably the M = 4 SYM CFT computation 
of these production rates suggested they are not affected by a hard IR-cut-off even 
for temperatures infinitesimally above the cut-off [11]. Intuitively this seems rather 
strange. At energies and temperatures close the QCD scale IR effects should start 
to affect the production rate. We shall find that for smoothly IR-cutoff AdS/QCD 
this is indeed the case. The robustness of our phenomenological result of how photon 
production rates are effected by changing the IR-cutoff is confirmed by a calculation 
by Mateos and Patino [12] of the photon production rate in AdS dual of a Af = 2 
theory with massive flavor. Here the flavor sector acts as the effective IR-cutoff, and 
we will be able to show this by relating the mass-parameter to the Soft- Wall cut-off 
scale. Soft-wall AdS/QCD is more crude than massive flavor models, of course, and 
this is evident in the lack of spectral peaks that we shall find. 

Photon production in a medium such as sQGP was discussed in detail both from 
strong and weak coupling point view in [11]. We briefly review this in section 2 and 
show there how the strong coupling calculation is modified by considering AdS/QCD 
instead of pure M = 4 SYM. In section 3, we present our solution and discuss its 
results in section 4 with a comparison to photon production in AdS duals of M = 2 
massive flavor theories. 



2. Photon and dilepton production 

One of the observational phenomena in RHIC is the spontaneous production of pho- 
tons from the sQGP of hot charged particles. This direct photon spectrum ought to 
be a good probe of the strongly coupled quark-gluon soup, as the weakly interacting 
photons should escape nearly unaffected from the small finite size collision area [13]. 

As is described in [11], we can therefore regard the dynamically formed sQGP 
to first approximation as a field theory at finite temperature. For a standard pertur- 
bative electromagnetic current coupling eJ™A M , the first order photon production 
rate is then given by [11, 14] 



1 The prediction of the critical deconfinemcnt temperature in soft-wall model is really quite 
surprising; The neglection of backreaction invalidates the thermodynamic relations that are used 
to derived it. Detailed studies that include the effects of backreaction are e.g. [22, 23, 24] for the 
backgrounds which are believed to correspond to large N c pure Yang-Mills theory and [25, 26] for 
the backgrounds that mimic QCD equation of state. 
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Here K = (k°, k) is a momentum 4-vector, ne(/c°) = l/(e^ k ° — 1) the Bose-Einstein 
distribution function, and the spectral density Xnv{K) is proportional to the imagi- 
nary part of the (finite temperature) retarded current-current correlation function 

X^{K) = -2 lm(G^(K)), 

G*f(K) = J d A Xe^ K - x (J^(0)Jr(X)),9(-x°) . (2.2) 

At finite temperature, Lorentz invariance is broken by the heat bath. We can use 
the remaining rotational symmetry plus gauge invariance to simplify the retarded 
correlator to 

G ^\K) = P^(K)U T (K) + P^(K)Tl L (K), (2.3) 

Here the transverse and longitudinal projectors are Pq (K) = 0, Pq^K) = 0, P?(K) = 
Sij - kikj/\k\ 2 , and P^,{K) = P^K) - P^ U {K), with i, j = x,y,z. We can trivially 
consider charged lepton production as well by considering non-lightlike momenta 
for off-shell photons: The leptons then result from virtual photon decay. Lepton 
pair production for each lepton species in the leading order of the electromagnetic 
couplings e and q, is given by [11, 14] 

dTfl = WfM^ [ ~ R2 ~ 4m2]V2( -^ 2 + ^)n b (k^(K)e(ky(-K 2 - 4m 2 ), 

(2.4) 

with e\ the electric charge of the lepton, m the lepton mass, 6(x) a unit step function, 
and the spectral density Xnv{K) is evaluated at the timelike momentum of the emitted 
particle pair. Note that both n T and U L contribute to the dilepton rate, but only n T 
contributes to the photon emission rate, because the longitudinal part must vanish 
for lightlike momenta, i.e. the unphysical longitudinal mode is not a propagating 
degree of freedom. 

Finally, fluctuation-dissipation relates the zero-frequency limit of the spectral 
density to the electrical conductivity a: 

o2 



a = lim ^-n B {W X>lu {k^ k = 0), (2.5) 



or, if k^ is lightlike 



e 2 



° = So W X,AK)\ rk]=k0 . (2.6) 



2.1 Photon and dilepton rates at strong coupling 

The AdS/CFT dictionary gives that the large N c limit of strongly coupled d = 
4 J\f = 4 SYM theory at finite temperature T has a dual description in terms of five 
dimensional AdS-supergravity in the background of a black hole [6] 

ds 2 = \-f( u )dt 2 + dx 2 + dy 2 + dz 2 ] + - f' , du 2 . (2.7) 
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Here f(u) = 1 — u 2 , with u G [0, 1] a dimensionless radial AdS coordinate related 
through u = (itTz) 2 to standard AdS coordinates, and R is the curvature radius of 
the AdS space. 2 The metric (2.7) has a horizon at u — 1 with Hawking temperature 
T and a boundary at u — 0. 

Qualitatively the same is expected hold for other 4-dim field theories. As a model 
for low energy QCD we shall take the AdS dual of chiral perturbation theory. This 
AdS/QCD consists of the fields A a L ^, A a Rfl , dual to the SU(N f ) L x SU(N f ) R currents 
and a scalar X dual to the quark condensate in an AdS background which is cutoff 
at some finite distance u = u [8]. To this we add an extra U(l) field, dual to the 
electromagnetic current J™. 3 Recall that Uq corresponds to the introduction of the 
QCD-scale in the field theory: it enforces the mass-gap by hand by explicitly cutting- 
off any dynamics in the IR. For the reasons we explained in the introduction, here 
we are going to use a soft wall cut-off [9, 10]. Formally we can introduce this cut-off 
by modifying the AdS bulk action to (we give only the term relevant for calculating 
the photon production rate) 



S ~ J Sx^g (- l -F AB F AB + ■ ■ =► S~ J d^e-* (-\f ab F ab + •••). 

(2.8) 

Here A,B = t, x, y, z, u and the "dilaton" takes the fixed form $ = cu where c = 

A 2 

, jffs , with Am the IR scale below which physics is cut-off. This introduction into 
the action is formal in the sense that (1) we shall not consider $ a dynamical field 
and (2) we assume that the presence of the cut-off does not affect the geometric 
AdS background, see also [10]. We thus still work with the metric (2.7) for the 
finite temperature version of AdS/QCD, but with the equation of motion for the 
fluctuations derived from action (2.8). We will discuss the validity of this approach 
in detail in section 4. 

For photon production, we need only the U(l) gauge field equation of motion 
9a (y/ge- cu g AB g CD F BD ) = with F AB = d A V B - d B V A the Maxwell field strength. 
The 4d electric fields are = F ti with % — x,y, z. Note that we use A as a vec- 
tor index and V B for the AdS gauge field. To compute the AdS boundary 2-point 
correlation function from which to extract the spectral density Xtw> we follow [11] 
and split the equation of motion into parts perpendicular (V x , V y = V±) and parallel 
{V z = V\\) to a predefined spatial three-momentum k = (0, 0, k), the Gauss constraint 
(Vq e.o.m.) and the radial AdS (V u ) equation of motion. After a Fourier transforma- 



2 We will keep to Lorentzian signature throughout since we seek information regarding the re- 
sponse of the thermal ensemble to small perturbations. This requires the use of real-time Green's 
functions [15]. 

3 In chiral pertubation theory this U (1) external field is a subgroup of the local SU (3) l x SU (3)_r 
[16]. 
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tion along t,x,y,z, and denning to = ^-f, q = o^f> we find respectively 



a?y± + 



f 



- c d u V ± + 



^ 2 - q 2 f 
uP 



V ± = 0, 



^-{qVt + u>V\\) - {d 2 u V t + i{2itT)ud u V u ) + c (d u V t + i{2itT)uV u ) = 0, 
uj 



UJ 



up 



{qV t +UjVn) + 



f 



- c)d u V\\ + d 2 u V\\ 



- i(2nT)q 



■d u f 



-c)V u + d u V u 



The equation of motion for V u , 

^ e -CU g UU (gtt dtFtu + ^HII 9|| F|| U ) = 0, 

can be simplified to 



V u 



(uduV t + qfduV l{ ) 



(2.9) 
(2.10) 

= 0. 
(2.11) 

(2.12) 
(2.13) 



2vrT (w 2 - g 2 /) 

Let us define E± = uuV± and E\\ = qVt + wV[|. From Eq. (2.9) and combining eq. 
(2.11) with eq. (2.10) and eq. (2.13) in the gauge V u = we obtain the two decoupled 
equations 



dlE, + 



f 

u 2 d u f 
P^ - g 2 /) 



c d u E ± + 



q 2 f 



uP 



E ± = 0, 



d u E\\ + 



q 2 f 



uP 



En = 0. 



(2.14) 
(2.15) 



We shall need to solve these two equations to obtain the spectral density x^y 

These 

differential equations (2.14) and (2.15) have three regular singular points at u — 
±1,0, and one irregular singular point at oo. 4 

Formal solutions for such equations are difficult to construct. Note that the 
irregular nature of the point at infinity becomes regular when we remove the IR- 
cutoff c. The irregular point, however, is outside the physical region of interest 
u E (0,1) and we can, for instance, solve the equations (2.9) and (2.15) near the 
boundary u — > using Frobenius expansion E = u x Xl^Lo a n.u n where the indicial 
equation has solutions for A = 0, 1. 

To solve the equations (2.14), (2.15) explicitly shall be the main part of this note. 
The solutions to these 5-d AdS equations of motion then give the 4-d field theory two 



4 Recall that an irregular singular point for a differential equation y" + P(x)y' + Q(x)y = is 
a point xq for which either \im x ^ Xa (x — xq)P(x) or lim x _ 5 . 2:o (x — xq) 2 Q(x) diverges. The point at 
infinity is irregular if lim a; ^ 00 (2 — xP(xj) or lim^oo x 2 Q(x) diverges. Using that / = (1 — u 2 ) 
one clearly sees how the introduction of c introduces a divergence in lim^ 
lim^oo 2 + 2i*7(l - u 2 ) + uc. 



, 2 - u(d u lnf - c) 
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point correlation as the functional derivative with respect to the boundary values of 
the on-shell AdS action 



S = j d 4 xdu ^e~ cu F AB F AB 



on-shell 



(2.16) 



with g 2 B = 16it 2 R/N 2 . Considering V u — gauge, we can write this as 

u=l 

\ - / ■ _ 

c 

°on-shell 



32ir 2 R 

N 2 T 2 /-oc 



d 4 x {^ge- cu V^) 



u=0 



16 



d 4 x e~ cu (V t d u V t -fVtduVi) 



u=l 



u=0 



(2.17) 



Fourier transforming to momentum space and selecting the particular direction cho- 
sen previously, we can rewrite the action using Minkowskian prescription formulated 
by Son and Starinets [15]. Together with the boundary condition that the solution 
of equations (2.9) and (2.15) must satisfy the incoming- wave boundary condition at 
the horizon u — 1, the resulting on-shell action becomes 



N 2 T 2 

^on-shell = hm 



doj dq 

16 uTo J (2tt) 2 ' 
/ 



./ 



q 2 f - u : 



d u E\\(u,K)E\\(u,-K) - ^d u E ± (u,K)E ± (u,-K) 



(2.18) 



From Eq. (2.18) and the condition described above, we can now compute the retarded 
current-current correlation function in term of two independent scalar functions 5 



K ' 8 u^o E\\(u,K) 

K J 8 u^o E ± {u,K) 



(2.19) 
(2.20) 



These functions in turn give us the photon and dilepton production at strong coupling 
via eq. (2.3) and eqs. (2.1) and (2.4). 



3. Solving the system 

In this section we will solve the equations (2.14) and (2.15) in order to compute 
the two scalar functions (2.19) and (2.20). Furthermore, we will take the imaginary 
part of those scalar functions and obtain the spectral density function (2.2) for finite 
temperature system. 

5 For a more detailed derivation of these functions see [17]. 
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The solutions which satisfy the incoming- wave boundary condition can be written 
in general as a Frobenius expansion near u — > 1 

E t (u) = (l-u)'^ 2 y t (u), (3.1) 

with yi{u) regular at u — 1. We will solve and discuss these equations extensively for 
lightlike momenta relevant for photon-production, both semi-analytically for asymp- 
totically small and large frequency and numerically for various values of the cut-off c 
for the full range of momenta. For timelike and spacelike momenta we only present 
the numerical solution. 



3.1 Lightlike momenta 

As has been explained in section 2, the longitudinal part of the scalar functions 
vanishes for lightlike momenta and we just need to compute the transverse part. 

3.1.1 Analytic solutions for lightlike momenta at low and high frequency 

We are mainly interested in the effect of the IR-cut-off on photon production as 
compared to the previous AdS photon production calculation for scale-invariant M = 
4 SYM [11]. In the low-frequency limit where its effect should be largest, we can solve 
(2.9) perturbatively using w C 1 as a small parameter. As noted in [11], there is a 
shortcut to do so. Given the two independent solutions 4>i ± ^02 to the differential 
equation 0" + A(x)cj)' + B(x)4> = 0, the Wronskian times exp(J x A(x')) is strictly 
conserved 

dsfef'AWlfat-Wj]) = 0. (3.2) 
The transverse scalar can be rewritten as 



IT (K) = lim IT (u, K) , 

u — ^0 



Il T (u,K) 



8 



e- cu {l-u z ) ' ' d u 



E ± (0,K) a E ± (0,K) 



(3.3) 



The imaginary part of the transverse scalar n T (w, K) is then propotional to the 
conserved Wronskian and therefore independent of the radial coordinate u: 

d u lm[U T (u,K)} = . (3.4) 

With this fact, we can evaluate the imaginary part of (3.3) at any given value of 
u which is convenient to our calculation. Let us choose u — 1. Because the trans- 
verse scalar (3.3) contains an explicit factor of (1 — u), only the pole in E±d u E± 
will contribute. Recalling that for any finite frequency oj the boundary conditions 
determine E±(u) to be of the form (3.1), we immediately see that the undetermined 
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regular part y contains no pole by definition. Therefore without needing to solve the 
equation motion we see that 

2" 



c 

8 



2e~ 



0(1) 



y(o) 



(3.5) 



The leading term in the limit uj <C 1 is the w-independent contribution to |y(l)/y(0)|. 
The determining equation (2.14) simplifies in that limit to effectively the first order 
equation (recall that u = q for lightlike momenta) 

d u d u E ± + (d u (\nf- cu))d u E = + O(co 2 ) . (3.6) 

The incoming wave boundary condition demands that the u = solution be regular 
at u — 1. Since / = (l — u)(l + u), this solution is the trivial constant one. Therefore 

n T (i, k) = !^L e - + o(u 2 ) . (3.7) 

8 

In Appendix A we compute the same answer directly by solving the differential 
equation perturbatively in u, which shows explicitly that E±(u) = constant + 0(uj) 
is indeed the correct solution to the boundary conditions. 

Given 11(1, K), the trace of spectral density function at low-frequency limit for 
lightlike momenta in photon production is proportional to its w-independent imagi- 
nary part 

X»(w = q) = -4 Im(n T (w = ,)) 

= ^+f(4 (3.8) 

For c = 0, we reproduce back the result from [11] at the first order. The vanishing 
of c corresponds to either the limit T — > oo or to removing the IR scale Air. We see 
explicitly our intuition confirmed that the trace of spectral density at low-frequency 
depends on the cutoff parameter c, while simultaneously reproducing the Af = 4 
result at high T. 

At high-frequencies we do not expect the IR-cut-off to have a major effect. Let 
us show that to leading order the spectral function is in fact independent of the value 
of c as one would expect. In this limit uj ^> 1, the argument leading up to eq. (3.5) 
does not hold 6 and one cannot obtain the answer without solving the equation of 
motion (2.14). Following [11], we will use the Langer-Olver method [18, 19] to find 
the solution. The first step is to redefine 

cu/2 

E ± (u) = —==y(u) (3.9) 

v-f( u ) 



6 Note e.g. that in the singular term (1 — u) %U) I 2 the order of limits u — > 1 and ui — >• oo do not 
commute. 
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for equation (2.14) and rewrite it as 

y"(x) = [u 2 H(x) + G(x)]y(x), 



(3.10) 



where H(x) = and G(x) = ^ - ^ - with x = -u e [-1,0]. For large 
oj the first term on the RHS dominates. Since it has a simple zero at x — 0, we can 
transform Eq. (3.10) to Airy's equation plus terms subleading in u. To do so, we 
introduce a new independent variable ( and change variables to 



c 



\dx 



H(x) 



x 



[i-x 2 y ' 

Choosing conditions ((0) = and ('(0) > determines ( to be 



(3.11) 



c = 



y/H{t)dt 



2/3 



(3.12) 



Rescaling y(x) to 



dx 



-1/2 



w , 



eq. (3.10) becomes 



d 2 W 



with 



5 + [AH(x)H"(x] - 5H' 2 {x)] + (G(x) 



(3.13) 



(3.14) 



(3.15) 



16C 2 lQH 3 (x) " H(x) 

For large u we may ignore V'(C) an d the equation reduces to Airy's equation. To 
leading order the solution is thus 



W(() = A Ai(cu 2 / 3 () + B Bi(u 2 / 3 () + ..., 



(3.16) 



The incoming-wave boundary conditions at the horizon imply that B should vanish. 
Thus the solution for E±(u) in asymptotic expansion for large u is 



E±(u) 



A e cu / 2 



— u 



1-1/4 



Ai(a; 2 / 3 C(-«)) + . . . , 



(3.17) 



and the transverse scalar at high-frequency limit equals 



N 2 T 2 



lim ( - H — d„ In 



u 



^Ai(a; 2 / 3 C(-^)) 
Ai(w 2 / 3 C(-«)) 



+ 



Before we move on, it is helpful to expand ((—u) around u = 



C(-«) = 



n2/3 



-i) 2 /V + 0( M 5 ) . 



(3.18) 



(3.19) 
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Therefore the middle term in (3.18), 



vanishes as u — > 0. Knowing the asymptotics of the Airy function the last term of 
(3.18) can be written as 

Ai(w 2 / 3 C(-w))) v > Ai(0) 

- 



and thus we obtain 



T _ N?T* (c e 2 -/a a; 2/33i/3 r(2/ 3 ) x 

11 " ~8~V2 + fa73) 

Note that this transverse scalar therefore depends on c. However, only the real part 
does. The trace of the spectral density function in high-frequency limit for lightlike 
momenta 



X » = -4 Im(IT r ) 

AT c 2 T 2 w 2 / 3 3 5 / 6 r(2/3) 



(3.23) 



4 T(l/3) 

does not depend on the cutoff parameter c at least up to first order and yields the 
same result as the calculation in J\f = 4 SYM. The fact that c does appear in the 
real part of the transverse scalar indicates that at first subleading order the spectral 
density function will likely differ from the M = 4 result. The numerical results in 
the next section bear this out. 

3.1.2 Numerical solution for lightlike momenta 

The analytic asymptotic solutions are a guidance to the full spectral function. The 
full solutions of equation (2.9) for non-zero c are very difficult to find, as we remarked 
earlier. This is due to the irregular singular point at u = oo for c ^ where analytic 
solutions are not known. In this subsection we are going to look for numerical 
solutions for non-zero c. 

We start from the general solution (3.1) which satisfies the incoming wave bound- 
ary condition. To set a parametrization of the initial conditions for the u — 1 regular 
function yi{u) = Ei(l — uf^l 2 of Eq. (3.1), we write the general solution as a poly- 
nomial expansion around u — 1, y(u) = Yl^=o an (^ ~ u ) n - Substituting (3.1) into 
equation (2.14) for lightlike momenta, we obtain the equation 

2 _ / 

\n—l 



^2 a n (n-i^j (1 - u) n 2 + c a n (n - i^j (1 - u) 1 

n=0 L 
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E 

m=0 



Om+1 



n-i- + — y — '- I ( I 



u) 



n+m— 1 



2m+2 



n+m— 2 



(3.Q4) 



The second sum (over m) arises from expanding = Y^=o 2^tt(1 — M ) n an d nj^p = 
^r!r(l — u) n . In order to find the coefficients a n , we have to solve this equation 
for each power of (1 — u) and obtain 

(1 — u)~ 2 : ao (arbitrary), 

! ^(c-1/2) 
2 1 — lui) 



(1-u) : a k = f k (u,c)a , 



(3.25) 



with /fc are functions of a; and c which vanish at cu = 0. This gives us y(u) and 
at u — 1 in terms of the above coefficients 



2/(1) = «o, 
y'(l) = 



-a - 



iw(c- 1/2) 
2(1 -iw) 



(3.26) 



These will be the two initial conditions for the differential equation for y(u). The 
explicit differential equation it must satisfy is 



y" + 



+ 



10O 



2u 



1 — u 1 — u 2 

2iu — u 2 



c y 



0J 2 U 



+ 



(1-u 2 ) 2 4(1 -u) 2 2(1 -u) \1 



2u 



+ c 



y = 0. 



(3.27) 



Notice that the initial conditions for y(u) still depend on an arbitrary constant a^. 
Physical quantities, such as the spectral density function, depend on ratios of y(u) 
and its derivatives and are independent of this constant. We are therefore free to set 
it to any value; we will choose a — 1. 

Let us express the trace of spectral density function in terms of y(u): 



N 2 T 2 



2 v 2/(0) ; 



(3.28) 



Alternately we could use the modified Wronskian formulation for H T ((K), eq. (3.3) 
and evaluate it at u — 1. An equivalent expression for the trace of spectral density 
function in this limit becomes 



uN 2 T 2 e _ c jy(l 



y(o)p 



(3.29) 
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3.1.3 The spectral density for lightlike momenta 

Solving eq. (3.27) numerically with initial conditions (3.26), we find the spectral 
density function f° r lightlike momenta for various values of the IR-cut-off c. 7 The 
results are shown in Fig. 1 and we clearly see the dependency at low frequencies on 
the IR-cut-off. The behaviour at high-frequency on the other hand appears less and 
less sensitive. What is remarkable is the similarity between this soft-wall AdS/QCD 



IF 
0.8" 




0.2 



1 \ 

1 2 3 4 5 

CO 

Figure 1: Trace of the spectral function for lightlike momenta in units of ^N 2 T 2 , 
plotted as a function of frequency with uj = k /(2irT). The solid line (red) shows 
the exact result for c = and the dashed lines downward show numerical analysis for 
c = 0.2, 0.419035, 0.5, 0.6, 1.5, 2, 3. 

result, Fig 1, for the trace of the spectral function for light-like momenta and of 
Mateos and Patino for massive flavor deformations of the AdS dual of M = 2 theories, 
Fig. 3 in [12]. As we will discuss in section 4, this similarity can be explained by 
relating the two computations. Inherently this then partially validates the soft-wall 
AdS/QCD model. 

There is, however, one fundamental difference between the result here and the 
massive M = 2 computation. Both models are thermodynamically unstable for large 
IR-cut-off, signalling the transition back to the confining regime. In the M = 2 model 
this is clearly illustrated by the appearance of thermal resonances in the spectral 
function when formally evaluated beyond the critical cut-off. Figl. shows that in 
AdS/QCD these resonances remain absent beyond the critical value c > 0.419035 
[10]. The absence of thermal resonances was presaged by Huot et al. [11]. Realizing 
that their results for photon production in the AdS dual of pure A^ = 4 SYM are 

7 Numerical solutions were obtained using the NDSolve routine in Mathematica. 
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unaffected by a hard-wall IR-cut-off, they speculated that this would be generic. It 
was premised on the fact that in the hard-wall case, the IR-cut-off is always inside 
the horizon. Rough dimensional analysis illustrates that the soft-wall case is similar: 
at the transition the cut-off scale c~ l ~ 2.5 is beyond the horizon u — 1. However, 
a similar argument holds for the massive j\f = 2 AdS dual. As we discuss in section 
4, the real reason for the absence of thermal resonances is probably simply that an 
blunt soft- or hard- IR-cut-off is too crude to capture this information. 

3.2 Timelike and spacelike momenta 

For time and space-like momenta, both I1 T and I1 L can contribute to spectral density 
function x^{K). Also for these cases, the mode equations (2.9) and (2.15) cannot be 
solved analytically for arbitrary frequency(w) and wave vector (q), and we determine 
the spectral function numerically. 

3.2.1 Numerical solution for transverse scalar function 

Following the same procedure in numerical analysis for lightlike momenta above, we 
substitute the general solution (3.1) for transverse direction into (2.9) and obtain an 
equation for y(u) 

ioo 2u 
1 — u 1 — V? 



&+{ — -—, 2- c )y'± 



oj - o (l - u ) 2iu — u iu ( 2u \ n , n nn . 

+ n 1 2 , 2 J yi~^ dl 2 +C ) ^ = - 3 - 30 

u{± — u A y 4(1 — u) z 2(1 — u) \1 — u A J\ 
As in the lightlike case, to determine the initial conditions we expand y(u) = 
J2^=o a n(^ ~ u ) n around u — 1, with 

a (arbitrary), a ± = ^ - a , a k = f k (u, q, c)a , (3.31) 

where again f k are functions of u, q and c which vanish at oo = q = 0. Using the 
modified Wronskian extension the imaginary part of transverse scalar function is 
therefore given by 

with y(u) a solution to eq. (3.30) with initial conditions determined from Eq. (3.31). 

3.2.2 Numerical solution for longitudinal scalar function 

Substitute (3.1) into the equation of motion for the longitudinal direction (2.15), we 
obtain 



„ ( iu 2uuf i , 

V W + [T^l ~ (1 -u 2 ){u 2 -q 2 {l -u 2 )) " C ) V W + 
2iu — u 2 iou ( 2uoj 2 



2 



q 2 (l-u 2 ) 



u(l — u 



2\2 



+ 4(1 - u) 2 2(1 - u) V(l - u 2 )(co 2 - q 2 (l - u 2 )) + ( 



2/|| = 0. (3.33) 
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Expanding y\\(u) = J2^=o a «(l ~ u ) n around u — 1 gives us 

u 2 -q 2 + tu 

a (arbitrary) , a x = ^ _ — L a^ , a k = f k (u,q,c)a Q , (3.34) 

where again fj. are functions of u, q and c which vanish at u = q = 0. The imaginary 
part of the longitudinal scalar function is 

with (u) the solution to (3.33) with initial conditions determined from Eq. (3.34). 



3.2.3 The spectral density for time- and space-like momenta 

Following formula (2.3), we can now write the trace of spectral function for time- 
and space-like momenta as 



ujN 2 c T 2 



-Jv±(i) 



\y± 



2 i i T 

2+ l + 2^ lm 



'yj(Q)' 
Co) 



(3.36) 



The complete results for are plotted in Fig. 2 and Fig. 3 as a function of frequency 
for several values of the spatial momentum. As we increase the value for c, one 
clearly sees that at low momenta the function decreases compared to c = 0. 




Figure 2: Spectral function trace XaV 6 ^ i n 
units of N 2 T 2 /2, plotted as a function of 
uj. The solid lines describe c = 0.5 and the 
dashed lines for c = while different col- 
ors represent q = 0(red), q = 1 (green), and 
q = 1.5(blue). 



Figure 3: Spectral function trace i n 
units of N 2 T 2 /2, plotted as a function of u 
for q = and various values of c = (black), 
c = 0.2(red), c = 0.5(green), c = 0.6(bue), 
c = 1.5(yellow), c = 2(magenta), and c = 
3 (cyan). 
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3.3 Electrical conductivity 

With the spectral density in hand, it is now straightforward to compute the electrical 
conductivity a. Here, we will use Eq. (2.6) as we have an analytic expression of the 
spectral density for lightlike momenta. Substituting (3.8) into (2.6) yields 

e 2 x£{u = q) 

e 2 N 2 exp(-c)k°(l + Q(k )) 
" fco™8^ fc°/T(l + 0(k°)) 
N 2 T 

= e 2 1 ^exp(-c), (3.37) 

with e the electric charge. We again note the presence of the scaling factor e~ c 
which dampens the IR-properties, including charge diffusion, of the system. Note 
in particular that this IR-suppresion is also present in the charge susceptibility S = 
iV 2 T 2 c/8(e c — 1) and the more "universal" diffusion constant D = o"/e 2 S = (1 — 
e~ c )/27rTc (see Appendix B). Physically this makes sense, as a mass-gap should 
dampen any hydrodynamic behaviour and the general AdS/CFT computation for 
scale-dependent currents 



~>AdS 



f d 4 xdu ^r g ^—F AB F AB , (3.38) 

J 9 e S\ U ) 

demonstrates this explicitly [20] 8 



D~ o t — n / d «»Sf («)••■ ( 3 - 39 ) 

9eS\ U ~ L ) JO 

4. Conclusion: Soft wall cut-offs as an IR mass-gap. 

The essential new ingredient in Soft- wall AdS/QCD is the ad- hoc cut-off of the radial 
AdS-direction. It is intended to capture the dominant effects of the scale dependence 
of QCD [8, 10]. However, its ad- hoc introduction opens it to criticism; especially 
when interpreted as a dilaton-profile without taking into account back-reaction effects 
or the dilaton equation of motion (see the footnote in the introduction). On the 
other hand the succesful results of the model [9, 10], suggest that it does capture the 
essential IR behaviour correctly. 

The result for AdS/QCD photon production supports this further. As previously 
emphasized it closely resembles photon production due to quarks for M = 2 theories 
with massive flavor in the probe approximation Nf <ti N c [12]. These theories descend 



8 This suggests a trivial violation of the PSS shear- viscosity-bound by IR-suppressing hydrody- 
namic behaviour. As the derivation of the viscosity in [20] suggests, however, and the explicit 
computation in massive J\f = 2 models shows [21], this is not case. 
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from brane-constructions in string theory, and therefore have no ad hoc component 
to criticise. Recall that in these theories, the probe approximation means that one 
may consider the flavor group as a global symmetry. The U(l) theory with respect 
to which photons are defined is a subgroup of this group and the tunable quark 
mass — a free parameter in the brane construction — functions as the scale in these 
theories. On the other hand, because the matter and symmetry content is different 
from QCD, one could question how relevant massive M = 2 SQCD results are to 
reality. The observation we make now is that the resemblence between the trace of 
the spectral function in these M = 2 SQCD theories as a function of the quark 
mass m and the AdS/QCD spectral function as a function of the IR-cut-off c can be 
mathematically explained. Both therefore demonstrate again that AdS / CFT results 
are remarkable universal and robust across fundamentally different theories. This 
is therefore strong support for soft-wall AdS/QCD, despite its ad-hoc IR-cut-off, as 
well as massive M = 2 SQCD, despite its unrealistic matter content, as descriptions 
of QCD. 

To relate the massive Af = 2 SQCD result to AdS/QCD, we note that Mateos 
and Patino showed that in M = 2 SQCD the defining equation relevant for the trace 
of the spectral function for lightlike momenta can be deduced from an action 9 

dudxodx! [-P(u)(d V ± ) 2 + /P(u)(9iVl) 2 + Q(u)(d u V ± ) 2 ] , (4.1) 



where 



P(u) 



u 



Q(u) = f 



uf 

(l-^,o(«)) 3 



1 u 3 fy/g^ mfl (u),u) u 2 f 2 (l - V4,oH) 3 
u uf u 6 g(ip mfi (u),u) 



(4.2) 



Here / = f(u) = (1 — u 2 ) is the non-extremality function in the D3-brane metric 
(2.7). The function ip m fi{u) is the solution to the embedding equation of motion for 
the D7-flavor brane derived from the DBI-action 

duy/g{i> m {u),u) = J du^(l - ^Wl-if>* + 4u 2 M 2 , (4.3) 

i.e. g(tp(u),u) is the induced metric on the flavor brane. The u = boundary 
behavior of the solution ip m $ = ^u 1 / 2 + A-u 3 / 2 + ... is determined by the masses m 
and condensate expectation value (qq) ~ A of the quarks. For the massless theory 



9 We only consider the D3/D7 brane set-up of [12]. The gauge/gravity duality for the D4/D6 
brane set-up they also consider is not yet fully understood. 
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■0 m= o,o = and ^fg = u~ 3 . Thus to find the spectral function, one must first solve the 
differential equation for ip m {u) with the appropriate boundary conditions and then 
solve the differential equation for V± [12]. The first step correctly incorporates the 
backreaction of the modified IR-physics as opposed to the AdS/QCD ad- hoc cut-off. 

The massive case ip m ,o(u) ^ is therefore a step more involved than the massless 
case, unlike AdS/QCD where the scale is a mild modification c ^ of the defining dif- 
ferential equation (2.14). However, searching for a closer match, one quickly realizes 

— * 

that the massless equation (for lightlike momenta u = k), 

d 2 u V ± + d u (\nQ)d u V ± + k 2 (l - f)^V ± = 

=► d 2 u V ± + d u (\n (f))d u V ± + k 2 (l - f)^pV ± = , (4.4) 

is exactly the AdS/QCD equation (2.14) for c = and we are therefore lead to 
consider a change of variables for the massive case that resembles that of the massless 
case. Thus we define a new variable u such that 



u 



\fg{^ m fl{u),u)) _ l 



du j = d u—~ (4.5) 

uf uf 

with / = f(u). By construction the parameter P in the new variable is identical to 
the massless case and Q is seen to be a mild modification 

P(u) ' 



uf(u) 



m = m *^*£ . (4,) 

Note that the solution to the massive embedding equation of motion, ip m ,o 7^ 0, is 
implicit in the transformation (4.5). In this new variable, however, we see, that its 
specific form only mildly modifies the massless differential equation 

(<i/)- ] 



ln(/)+ln((l-< )^ 



u 



a a V ± + ^(l-f) j^^ V ± = 0(4.7) 



u(u) ' 

and the close relation to AdS/QCD is now apparent. The resemblance of the spectral 
functions is especially explained, if we recall that it is primarily determined by the 
u — behaviour of the solution (2. 20). 10 As we know what the u = behaviour 
of the solution ip mt0 = ^u 1 / 2 + . . . must be, Eq. (4.5) shows that asymptotically 
u = u + ^u 2 + . . . and we can putatively identify the mass m with the IR-cut-off c: 



~ 2 2 

cu ~ ln(l-^ i0 ) 3 -=ln(l-— « + ...) 3 -ln(l -_« + ...) 
5 

~ m 2 u + ... (4.8) 

4 



10 Onc should be careful in that the change of coordinates (4.5) in principle will also change the 
boundary conditions one must impose. 
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The map between AdS/QCD and M = 2 SQCD is not exact; clearly we should 
not have expected it to be. The latter shows thermal resonances in the spectral 
function for masses m > 1.3092 which is the value beyond which the AdS black- 
hole solution becomes thermodynamically unstable [12]. The AdS/QCD description 
is much cruder as is no resonances show up even beyond the unstable regime c > 
0.419035. These thermal resonances are encoded in the subtleties of the embedding 
function ip m ,o(u) which carries more information than just the mass as an IR-cut- 
off. Precisely, the embedding function determines whether the flavor D7-brane is 
in "Minkowski embedding" or "black hole embedding" corresponding to the low T 
confining or high T deconfining phase [12]. Clearly, the M = 2 SQCD theory has a 
more detailed description at the physics. On the other hand, the results here do show 
that in the stable phase the simple AdS/QCD model describes the IR-consequences 
of a mass-gap remarkably well and the above derivation explains mathematically 
why. This in itself lends support to continue to study AdS/QCD as a good toy 
model for real-world physics. 
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A. Spectral function low frequency limit for lightlike momenta 

Here we find an analytic expression for the low-frequency limit of the transverse 
scalar and spectral density for lightlike momenta by solving the differential equation 
for the E±(u) perturbatively, rather than using the Wronskian shortcut, explained 
above eq. (3.2). 

We first extract the other regular singularity at u — — 1, writing 



with Y(u) regular at u — 1 and substitute this into (2.9). Changing variables to 
v = 1/2(1 — u), we obtain the differential equation 



In the absence of the IR-cutoff, c = 0, we recognize a hypergeometric equation with 
solution [11] 



E ± (u) = (1 - u)~ iul/2 {l + u)-^ 2 Y{u), 



(A.l) 



v(l — v)Y" + [(1 — iu) — (2 — iu — u — 2c)v - 2cv 2 ] Y' 

— | - [—u — icu + ico 2 ~\ — c[uv — iu + ituv] 1 Y — 0. (A.2) 
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As we noted earlier, the presence of c changes the nature of the equation and no 
formal solution is known. On physical grounds we expect the effects of c to dominate 
the low frequency part of the spectral function. Expanding Y(u) as 

Y = Y + uY 1 + cu 2 Y 2 + u 3 Y 3 + ■■■, 



we find to first order in ou, 



v(l - v)Y£ + [l-2v + 2cv(l - v)\Y^ = 0, 



u 1 : v{\ - v)Y(' + [v - i(l - v)\Y^ + [1 - 2v + 2cv(l - v)]Y{ 



+ {|(1 + i) + c[v - i(l - v)]} Y = 0. 



These two equations have solutions 



Y (v) = A + B [e" 2c Ei(2c - 2cv) - Ei(-2cu)] , 



(A.4) 

(A.5) 
(A.6) 

(A.7) 



Y^v) = C + -pn(u-l)+iliiu] + 

[e" 2c Ei(2c - 2cv) - Ei(-2cu)] 



D + — [\n(v - 1) +i\nv] 



with A, B, C, D constants of integration and Ei(x) = — ^j-dt the exponential 
integral function. To determine the integration constants, recall that by construction 
the solutions must be regular as v — > (u — > 1). Since the exponential integral Ei(v) 
diverges at v — 0, we must set B = 0. To determine regularity of Yi(v), recall that 
Ei(x) can be written as 

(-l) n x n 



Ei(— x) = 7 + In a; + 



n=l 



nln 



for x > 0, 



(A.8) 



with 7 the Euler-Mascheroni constant. Since the variable v G [0,1/2], and c > 0, 
regularity at v — demands D = iA/2. For convenience, let us also redefine the 
constant C = iCA/2. Substituting those constants into Yi, we obtain the solution 
for E± in the low frequency limit 



E ± (u) = A(l - u)- w /\l + u)-^ 2 |l + i| \ C + e- 2c Ei(c(l + «)) - Ei(c(u - 1)) 



-i In 



tt + 1 



+ In 



1 



u 



(A.9) 



Using the definition of the exponential integral function, we straightforwardly obtain 
the leading low-frequency contribution to transverse scalar function 



N 2 T 2 

n> = q ) = - ! ^4- 



icuN 2 T 2 



16 



CO co too , ... „ , . 

~ % 2 ~ 2 + Y (ce Ei ( c ) - cEi (~ c ) 

2/ - (e" c + e" c ) + G(o; 2 )] . 



1) + 0{u 2 ) 



(A.10) 



This is the exact answer. The imaginary part computed via the conserved Wronskian 
shortcut (3.7) clearly agrees. 
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B. The susceptibility and the diffusion constant 



We follow the procedure to compute the diffusion constant described in [7]. Using 
the gauge V u = 0, we can rewrite equation (2.10) as 



Mi 



uf„„ uf q 
— V. — c — V, Vt- 



quo 



quo 



(B.l) 



Substituting into equation (2.13) we obtain a second order differential equation for 



E = y; 



E" + 



(uf)' 

uf 



E' 



'oo 2 - q 2 f (uf)' 



E = 0. 



(B.2) 



uf 2 uf 

Imposing the same incoming-wave boundary condition as before and extracting the 
singularity at the horizon u — 1, we rewrite E = (1 — u)~ luJ / 2 y, where y is a regular 
function at the horizon. The function y must obey the equation 



y" + 



iu (uf)' 
1 — u 



uf 



y' + 



+ 



tu(tu + 2) tu((uf)' - cuf) u 2 -g 2 f (uf)' 



4(1 -u) 2 2uf(l-u) 



up 



uf 



y = 0. 



(B.3) 



For low frequency and momentum, we again solve the equation pertubatively in u 
and q 

y(u) = yoo + uyw + q 2 yo2 h — • (B.4) 

Up to first order in oj and q 2 , we find the system of equations 



<A° : 2/oo + 



toY : y'( + 



cA 2 : y'02 + 



(uf)' 



uf 



i (uf)' 

yoo - c —ryoo = °> 



uf 



1 — u 



y'oo + 



(uf)' 
uf 



2/io + 



i((uf)' - cuf) 



(uf)' 
uf 



(uf)' 

y02 - C ^02 

uj 



2(1 -u) 2 2uf(l-u 

f 



yoo 



uf)' 
-c —J/10 = 0, 



uf 



up 



:2/oo = 0. 



(B.5) 



Using the same analysis for the low frequency of spectral function as described in 
the previous Appendix, the solutions regular at u — 1 are found to be 

2/oo = Ae cu 

j A 

y w = — e cu+c [Cio + 2Ei(-cu) - e c Ei(-c(l + u)) - e~ c (Ei(c(l - u)) - \n(u - 1))] 
A 



2/02 = — e cu+c [C 2 - 2Ei(-cu) + e c Ei(-c(l + u)) 



2c 



+e- c (Ei(c(l - «)) + 2 In u - ln(w 2 - 1))] , 



(B.6) 
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where A and C±o, C02 are constants independent of u. We can determine A in terms 
of the boundary values of V t and V\\ at u — >■ defined as 

limVi|(u) = V?. (B.7) 

Substituting the solution for E = V' t into equation (B.l) and taking limit u — > 0, the 
integration constants C w , C 02 drop out and we can determine A to be 

q 2 V t ° + wgVf 

A = itoe c - f (1 - e- c ) g 2 + 0(w 2 , wg 2 , g 4 ) ' (B ' 8) 

We recognize the hydrodynamic pole and as explained in [7] we can now compute 
the time-time component of the retarded thermal Green's function of two currents 

8{tuj - - > q 2 ) 

Thus the time-time component of the spectral density function at low frequency and 
momentum equals 

^,^-2^ = ^,1^ ,..., (B,0, 

(l-e- c ) 

with D = v 27tTc ' the diffusion constant. Comparing the result with the universal 
hydrodynamic behaviour 

x ^ rk) = ^^f + -' (B11) 

the charge susceptibility S is seen to equal S = ^ and naturally satisfies the 
Einstein relation S = a/e 2 D. 
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